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^H Abstract 

O 

We identify a class of supersymmetric SU{A)c x SU{2)l x SU{2)r models in which 

c3 imposing essentially perfect t-h-r Yukawa coupling unification at Mqut yields a mass 

^ close to 122-126 GeV for the lightest CP-even (SM-like) Higgs boson. The squark 

r^^ and gluino masses in these models exceed 3 TeV, but the stau and charginos in some 

CN cases can be considerably lighter. We display some benchmark points corresponding 

,__i to neutralino-stau and bino-wino coannihilations as well as A-resonance. The well- 

i-^ known MSSM parameter tan/3 is around 46-52. 
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1 Introduction 

In a recent paper, hereafter referred to as [1], we explored the LHC imphcations of a 
supersymmetric SO(IO) model with third family Yukawa coupling unification [2, 3, 4] 
and 5*0(10) compatible gaugino mass ratio at Mqut given by M3 : M2 : Mi = 
2 : -3 : -1 [5]. Here M3, M2 and Mi respectively stand for SU{3)c, SU{2)l and 
U{1)y soft supersymmetry breaking (SSB) gaugino masses. It was also assumed in 
[1] that the SSB mass^ terms for the up and down type Higgs doublets of the minimal 
supersymmetric standard model (MSSM) are universal at Mqut- It was shown in 
[1] that t-h-T Yukawa coupling unification and neutralino dark matter abundance are 
readily compatible with each other in this S'O(IO) model, and a prediction of 122 - 
124 GeV for the lightest CP-even Higgs mass (with a theoretical uncertainty of ±3 
GeV [6]) was highlighted. 

Spurred by these results, we investigate here t-h-r Yukawa coupling unification 
in the framework of supersymmetric S'?7(4)c x SU{2)l x SU{2)pi (4-2-2, for short). 
The 4-2-2 structure [7] allows us to consider non- universal gaugino masses, while 
preserving Yukawa coupling unification and universality of the up and down type 
SSB Higgs masses. Previously, t-b-r Yukawa unification in the framework of 4-2-2 
was investigated in [8, 9, 10], where non-universal up and down type SSB Higgs masses 
were assumed. An important conclusion reached in [8, 9] is that with same sign non- 
universal gaugino soft mass terms and /x > 0, where /i denotes the coefficient of the 
bilinear Higgs mixing term, Yukawa unification in 4-2-2 is compatible with neutralino 
dark matter, with gluino co- annihilation [8, 9, 11, 12] being the unique dark matter 
scenario. 

By considering opposite sign gauginos such that /i < 0, M2 < and M3 > 0, it 
is shown in [10] that Yukawa coupling unification consistent with the experimental 
constraints can be implemented in 4-2-2. Yukawa coupling unification in this case is 
achieved for mig > 300 GeV, as opposed to niiQ > 8 TeV for the case of same sign 
gauginos. The finite corrections to the b-quark mass play an important role here [10]. 
With M2 < 0, M3 > and /i < in 4-2-2, we can also obtain additional contributions 
with the correct sign for the muon anomalous magnetic moment [13]. This enables 
us to simultaneously satisfy the requirements of t-b-T Yukawa unification, neutralino 
dark matter abundance and constraints from {g — 2)^, as well as a variety of other 
bounds. 

A crucial difference between this paper and previous studies of 4-2-2 with t-b-T 
Yukawa unification is the universality of the up and down type SSB Higgs masses 
assumed here. As we will see, the condition mj^^ = m\j at Mqut significantly reduces 
the range of allowed masses for the lightest CP-even Higgs boson. With rn^j 7^ rn^j 
and imposing t-b-T Yukawa unification in 4-2-2, one finds 80 GeV < nih < 128 GeV. 
However, in our present model, this range shrinks to 122 GeV < mh < 126 GeV [10], 
which is in good agreement with the recently reported evidence from the ATLAS and 



CMS experiments [14, 15] for a Higgs boson of mass around 125 GeV. 

The outline for the rest of the paper is as follows. In Section 2 we briefly describe 
the SSB terms in the 4-2-2 model. In Section 3 we summarize the scanning procedure 
and the experimental constraints that we have employed. In Section 4 we present our 
results, focusing in particular on the mass of the lightest (SM-like) CP-even Higgs 
boson. The squarks and gluino turn out to be relatively heavy, > 3 TeV. However, the 
stau and chargino in some cases can be considerably lighter. The table in this section 
highlights some benchmark points which can be tested at the LHC Our conclusions 
are summarized in Section 5. 

2 Supersymmetric 4-2-2 and SSB Terms 

In 4-2-2 the 16-plet of SO{10) matter fields consists of ip (4, 2, 1) and ipc (4, 1, 2). The 
third family Yukawa coupling ifjcipH, where H(l,2,2) denotes the bi-doublet (1,2,2), 
yields the following relation at Mqut [2]: 

Yt = Y, = Y, = Y,^. (1) 

Supplementing 4-2-2 with a discrete left-right (LR) symmetry [7, 16] (more precisely 
C-parity) [17] reduces the number of independent gauge couplings in 4-2-2 from three 
to two. This is because C-parity imposes the gauge coupling unification condition, 
9l = Qr at Mqut- We will assume that due to C-parity the SSB mass terms, induced 
at Mqut through gravity mediated supersymmetry breaking [18], are equal in mag- 
nitude for the squarks and sleptons of the three families. The tree level asymptotic 
SSB gaugino masses, on the other hand, can be non-universal from the following 
consideration. From C-parity, we can expect that the gaugino masses at Mqut as- 
sociated with SU{2)l and SU{2)r are the same (M2 = M!^ = M^). However, the 
asymptotic S'f/(4)c and consequently SU{3)c SSB gaugino masses can be different. 
With the hypercharge generator in 4-2-2 given by F = a/2/5 {B — L) + a/3/5 I3R, 
where B — L and I^r are the diagonal generators of SU{4:)c and SU{2)r, we have the 
following asymptotic relation among the three SSB gaugino masses: 

Ml = ^M2 + ^M3. (2) 

The supersymmetric 4-2-2 model with C-parity that we consider thus has two inde- 
pendent parameters (M2 and M3) in the gaugino sector. 

As was pointed out in ref [1], with non- universal gaugino masses at Mqut, ra- 
diative electroweak symmetry breaking (REWSB) compatible with perfect or near- 
perfect t-b-T Yukawa unification can occur even if we set rriH^ = f^H^- Note that this 
is quite difficult, if not impossible, to achieve if gaugino mass universality is imposed 
at Mgut [19]. 



The fundamental parameters of the 4-2-2 model that we consider are as follows: 



mi6, mio, M2, M3, v4o, tan/3, sign(/i). (3) 

Here rriiQ is the universal SSB mass for MSSM sfermions, rriio is the universal SSB 
mass for MSSM Higgs doublets, Aq is the universal coupling of the SSB terms for the 
trilinear scalar interactions (with the corresponding Yukawa coupling factored out), 
tan/3 is the ratio of the vacuum expectation values (VEVs) of the two MSSM Higgs 
doublets, and the magnitude of /i, but not its sign, is determined by the REWSB 
condition. Although not necessary, we will assume that the gauge coupling unification 
condition Qs = Qi = g2 holds at Mqut in 4-2-2. Such a scenario can arise, for example, 
from a higher dimensional 5*0(10) model [20], or from a SU{8) model after suitable 
compactification [21]. 

3 Experimental Constraints, Scanning Procedure 
and Parameter Space 

We employ the IS A JET 7.80 package [22] to perform random scans over the funda- 
mental parameter space. In this package, the weak scale values of gauge and third 
generation Yukawa couplings are evolved to Mgut via the MSSM renormalization 
group equations (RGEs) in the DR regularization scheme. We do not strictly enforce 
the unification condition gs = gi = (72 at Mqut, since a few percent deviation from 
unification can be assigned to unknown GUT-scale threshold corrections [23]. The 
deviation between gi = g2 and g^ at Mqut is no worse than 3 — 4%. For simplicity, we 
do not include the Dirac neutrino Yukawa coupling in the RGEs, whose contribution 
is expected to be small. 

The various boundary conditions are imposed at Mqut and all the SSB parame- 
ters, along with the gauge and Yukawa couplings, are evolved back to the weak scale 
Ml- In the evolution of Yukawa couplings the SUSY threshold corrections [24] are 
taken into account at the common scale Msusy = y/i^JW^^ where mi^ and mi^ de- 
note the masses of the third generation left and right-handed stop quarks. The entire 
parameter set is iteratively run between Mz and Mqut using the full 2-loop RGEs 
until a stable solution is obtained. To better account for leading-log corrections, 
one-loop step-beta functions are adopted for the gauge and Yukawa couplings, and 
the SSB parameters rui are extracted from RGEs at multiple scales rrii = mi{mi). 
The RGE- improved 1-loop effective potential is minimized at Msusy, which effec- 
tively accounts for the leading 2-loop corrections. Full 1-loop radiative corrections 
are incorporated for all sparticle masses. 



The approximate error of ±3 GeV in the ISAJET estimation of the hghtest CP- 
even Higgs boson mass largely arises from theoretical uncertainties [6] in the calcu- 
lation and to a lesser extent from experimental uncertainties. 

An important constraint on the parameter space arises from limits on the cos- 
mological abundance of stable charged particles [25]. This excludes regions in the 
parameter space where charged SUSY particles become the lightest supersymmetric 
particle (LSP). We accept only those solutions for which one of the neutralinos is the 
LSP and saturates the WMAP bound on relic dark matter abundance. 

We have performed random scans for the following parameter range: 

< mi6 < 6 TeV 

< mio < 5 TeV 

< M3 < 5 TeV 

-5TeV < M2 < 

35 < tan /3 < 55 

—3 < Ao/miQ < 3 

/i < (4) 

We set rrit = 173.3 GeV [26], and we show that our results are not too sensitive to one 
or two sigma variation from the central value of rrit [27]. Note that mb{mz) = 2.83 
GeV, which is hard-coded into ISAJET. 

In order to obtain the correct sign for the desired contribution to {g — 2)^, we will 
focus here on the case /x < and M2 < 0. 

In scanning the parameter space, we employ the Metropolis-Hastings algorithm 
as described in [28]. The data points collected all satisfy the requirement of REWSB, 
with the neutralino in each case being the LSP. After collecting the data, we impose 
the mass bounds on all the particles [25] and use the IsaTools package [29] to imple- 
ment the various phenomenological constraints. We successively apply the following 
experimental constraints on the data that we acquire from ISAJET: 

rrih (lightest Higgs mass) > 114.4 GeV [30] 

BR{Bs -^ fi+fi ) < 4.5 X 10-9 [31] 

2.85 X 10-^ < BR{b -^ S7) < 4.24 x 10"^ {2a) 

0.15 < ^S^"^-^^^"^ < 2.41 (3a) 

f^CDM/i^ = 0.1123 ± 0.0035 (5a) 
< A(^ - 2)^/2 < 55.6 X 10"^° [13] 
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Figure 1: Plots in R — mig, -R — tan (3, R — Aq/uiiq and R — uiiq planes. Gray points 
are consistent with REWSB and neutralino LSP. Green points satisfy particle mass 
bounds and constraints from BR{Bs — ?■ fi'^fi^), BR{b -^ 57) and BR{Bu — ?■ tz/,-). In 
addition, we require that green points do no worse than the SM in terms of {g — 2)^. 
Brown points belong to a subset of green points and satisfy the WMAP bounds on 
neutralino dark matter abundance. 

4 Sparticle Spectroscopy 

In order to quantify Yukawa coupling unification, we define the quantity R as. 



R = 



min(yj,|/b,l/r) 



(5) 



In Figure 1 we show the results in the R — rriiQ, R — Ao/miQ and R — niio planes. 
The gray points are consistent with REWSB and neutralino LSP. The green points 
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Figure 2: Plots in M2 — M1, yU — Mi, M3 — rriie and miQ — miQ planes. The color coding 
is the same as in Figure 1. In addition, we have used yellow color to denote a subset of 
the green points, that have Yukawa unification better than 5%. Brown points belong 
to a subset of yellow points and satisfy the WMAP bounds on neutralino dark matter 
abundance. 



satisfy particle mass bounds and constraints from BR{Bs — )■ ji^ji^), BR{h — )■ 57) 
and BR{Bu — ?■ tv^.). In addition, the green points do no worse than the SM in terms 
of [g — 2)^. The brown points belong to a subset of the green points and satisfy the 
WMAP bound on neutralino dark matter abundance. 

In the R — rriiQ plane of Figure 1, we can see realization of perfect Yukawa unifica- 
tion consistent with all constraints mentioned in Section 3. This is possible because 
we can implement Yukawa unification for relatively small values of rriiQ (~ ITeV). 
A similar conclusion was reached in ref. [10] but with relatively lighter values for 
rriiQ. This difference is related to the additional condition rriH^ = fT^H^ which we have 
imposed in the current paper. 
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Figure 3: Plots in the R — rrih plane. The color coding is the same as in Fig 1. In 
the left panel, the top mass is set equal to its current central value (m^ = 173.3 GeV 
[26]), while la deviation in the top mass {mt = 17AA GeV) is considered in the right 
panel to investigate the effect on the CP-even light higgs mass. 



Note that there is more than an order of magnitude reduction in the rriiQ values 
required for Yukawa unification, as compared with the case fi > and universal 
gaugino masses [4, 8]. There is upper bound for rriiQ once Yukawa unification is 
imposed. Indeed, one predicts tan (3 ^ 47 for perfect Yukawa unification. In i? — 
^o/'^ie plane, we see that perfect Yukawa unification can occur for Aq < 0, in 
contrast to the SO(IO) case with universal gaugino mass condition where perfect 
Yukawa unification takes place for Aq/ttiiq ~ 2.6. From the R — rriiQ plane we see 
that perfect Yukawa unification also puts an upper bound on rriiQ. 

In Figure 2 we show the results in the M2 — Mi, /i — Mi, M3 — rriiQ and miQ — rriiQ 
planes. The color coding is the same as in Figure 1. In addition, we have used yellow 
color for the subset of green points satisfying Yukawa unification to better than 5%. 
The brown points belong to a subset of green points which satisfy the WMAP bounds 
on neutralino dark matter abundance. 

From the M2 — Mi plane, we see that some of the brown points lie near the line 
corresponding to M1/M2 ~ — 2, which yields at low sale M^ ~ M^, the condition 
for bino-wino coannihilation [34] scenario. At the same time there is no colored point 
near the line corresponding to M1/M2 ~ 2, which means that in our model we do not 
have bino-wino mixed dark matter [35]. 

In the /i — Ml plane we have shown lines corresponding to fi/Mi = 1 and /u/Mi = 
— 1. The yellow points describing 5% or better Yukawa unification are far from these 
lines, which means that we do not have bino-higgsino mixed dark matter in this 
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Figure 4: Plots in the rrif — rrih, itia — tTih, '>Tij^± — rrih, mi — rrih, mi — m^ and 
iTLg — rufi planes. The color coding is the same as in Figure 1. The brown points in 
the TUf — TUh plane represent stau-neutralino coannihilation, A-funnel region in the 
m^A — rrih panel, bino-wino dark mater and bino-wino coannihilation in the m^± — m^h 
plane. In the m.i — m,h, m,i — m,h and m,g — m^h planes, the blue points form a subset of 
the green points and do not satisfy Yukawa unification. Blue points represent stop- 
neutralino coannihilation in the m^ — m^h plane, and gluino-neutralino coannihilation 
in the m,g — m,h plane. 



model. Accordingly, as shown in ref [36], it will be difficult to find dark matter 
predicted in this model in direct and indirect searches. 

The M3 — mi6 panel indicates that the minimal value of M3 which corresponds to 
5% or better Yukawa unification (yellow points) is above 1.5 TeV, which means that 
the colored particle are much heavier than 1 TeV. In particular, the minimum gluino 
mass is close to 3 TeV. 

From the rriiQ — rriiQ plane we see that one can reduce number of independent 
parameters, assuming mie = mio, and still retain decent Yukawa unification with the 
correct relic abundances. 

In Figure 3, we present results in the R — ruh plane for rut = 173.3 GeV and 
rrit = 174.4 GeV. The color coding is the same as the one used in Figure 1. In the left 
panel, the top quark mass is set equal to its central value {rrit = 173.3 GeV [26]), and 
it shows that we can find solutions with rrih ~ 125 GeV consistent with the constraints 
mentioned in Section 3 as well as t-b-r Yukawa unification. It is interesting to note 
that the Higgs mass interval is significantly reduced compared with the 4-2-2 model 
with rriH^ 7^ ?^//^ where it lies in the interval 80 GeV < rUh < 129 GeV, compatible 
with perfect Yukawa unification condition. In the current model, as seen from Figure 
3, the Higgs mass lies in the range 118 GeV < rrih < 127 GeV (gray points). This 
is further reduced once collider bounds and precise Yukawa unification are imposed 
(green points) and the range becomes 122 GeV < rrih < 126 GeV, which is in very 
good agreement with the recent ATLAS and CMS measurements [14, 15]. 

As far as neutralino dark matter abundance compatible with the WMAP bound 
is concerned, in Figure 4 we present results in the m^ — rrih, "^a — "^/i, "^x='= ~ "^/i' 
m^ — rrih, "^f, — "^/i and rrig — rrih planes. The color coding is the same as in Figure 1. 
The brown points in the mf — rrih plane represent only stau-neutralino coannihilation, 
they also represent A- funnel solution in the rriA — rrih panel, as well as bino-wino dark 
matter and bino-wino coannihilation in the m^± —rrih plane. In the mi — rrih, TT^i~^T^h 
and rrig-mh planes, the blue points are a subset of the green points and do not satisfy 
Yukawa unification. 

If the preliminary evidence for a Higgs with rrth = 125 GeV from ATLAS/CMS is 
verified, we see from the rrif — rrih panel solutions with 200 GeV < rrif < 600 GeV, 
linked to stau-neutralino coannihilation. Similarly, the rriA — rrih plane shows A- 
funnel solutions with 500 GeV < tua < 1200 GeV. The m^± — rrih panel shows 
that we can have bino-wino dark matter, and bino-wino coannihilation is realized in 
nature with 100 GeV < m^± < 600 GeV. We did not find solutions with neutralino- 
stop coannihilation compatible with 5% or better Yukawa unification. However, such 
solutions appear if the Yukawa unification condition is relaxed (blue points). In the 
"fT^b — fT^h plane we do not see either brown or blue points which shows that in the 
4-2-2 model it is not easy to find neutralino-sbottom coannihilation, a scenario which 
is easily realized in SU(5) [37]. We also did not find neutralino-gluino coannihilation 
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compatible with Yukawa unification (no brown points in the rrig — rrih plane), but 
there are blue points which show that we can have neutralino-gluino coannihilation 
for NLSP gluino up to 1.5 TeV mass, if we ignore Yukawa unification. 

Finally, in Table 1 we present three benchmark points which satisfy the various 
constraints mentioned in Section 3. Point 1 displays a solution corresponding to 
stau-neutralino coannihilation, point 2 depicts an A-funnel solution, point 3 stands 
for a solution with bino-wino coannihilation, and point 4 corresponds to a solution 
with rrtiQ = rriio. Each point in Table 1 has the CP-even lightest Higgs boson mass 
of about 125 GeV. 

5 Conclusion 

We have reconsidered t-b-r Yukawa unification in supersymmetric SU{A)c x SU{2)l x 
SU{2)fi within a slightly revised framework in this paper. The main difference from 
most previous investigations stems from the assumptions we make related to the soft 
supersymmetry breaking parameters. We set the masses of the two MSSM Higgs dou- 
blets to be equal at Mqut- The ramifications of these slightly different assumptions 
for TeV scale physics turn out to be quite startling, with the low energy predictions 
being quite different from previous studies. In particular, in the present framework 
with fi and M2 both negative, demanding perfect or near-perfect Yukawa unification 
yields a Higgs boson mass of 122 - 125 GeV (with an uncertainty of ±3 GeV). 

The solutions, obtained using the ISAJET software, are compatible with all ex- 
perimental observations, as well as the WMAP dark matter constraint. The masses of 
the gluino and first two family squarks are found to lie in the 2.7-5 TeV range, while 
the lightest stop (top squark) weighs at least 2 TeV or so. The MSSM parameter 
tan/3 is around 46 - 52. 
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Point 1 


Point 2 


Point 3 


Point 4 


mi6 


1729 


2777 


2406 


2016 


Ml 


-809 


1262 


1177 


-1046 


M2 


-4180 


-685 


-633 


-4867 


M3 


4247 


4183 


3892 


4686 


mio 


965 


640 


227 


2016 


tan/3 


49.2 


48.6 


48.2 


48.2 


Ao/mo 


1.04 


-2.74 


-2.82 


-2.81 


^t 


173.3 


173.3 


173.3 


173.3 


/i 


-3762 


-5765 


-5305 


-4539 


A(^ - 2), 


0.62 X 10-1° 


0.50 X 10-1° 


0.59 X 10-1° 


0.43 X 10-1° 


rrih 


125 


125 


125 


125 


rriH 


949 


1160 


1197 


1048 


rriA 


943 


1152 


1190 


1042 


mH± 


954 


1164 


1201 


1053 


A,l 2 


405, 3608 


556, 663 


515, 614 


508, 4188 


rriyo 


3769, 3795 


5739, 5740 


5283, 5283 


4540, 4552 


m-± 


3605, 3812 


666, 5738 


616, 5282 


4193, 4555 


m^ 


8496 


8455 


7891 


9340 


"^n^.fl 


7828, 7386 


7691, 7705 


7128, 7140 


8674, 8138 


"^il,2 


6265, 6943 


5978, 6209 


5579, 5800 


6541, 7387 


^d^,« 


7828, 7388 


7692, 7706 


7128, 7141 


8675, 8139 


^61,2 


6430, 6903 


6128, 6225 


5716, 5813 


6743, 7346 


m^ji 


3147 


2807 


2435 


3650 


m^3 


2909 


2344 


2037 


3376 


^eL,R 


3155, 1743 


2810, 2810 


2439, 2439 


3658, 2040 


"^ri,2 


426, 2925 


1617, 2330 


1388, 2028 


710, 3402 


f^s/(pb) 


0.67 X 10-11 


0.11 X 10-11 


0.11 X 10-11 


0.36 X 10-11 


o-sD(pb) 


0.46 X 10-9 


0.30 X 10-1° 


0.45 X 10-1° 


0.20 X 10-9 


^CDMh'^ 


0.26 


0.11 


0.11 


1.22 


R 


1.03 


1.02 


1.04 


1.03 



Table 1: Benchmark points for the 4-2-2 model with /i < 0. Point 1 shows stau- 
coannihilation, point 2 represents the A-resonance solution, point 3 depicts bino-wino 
coannihilation, and point 4 displays a solution with rriiQ = itliq. 
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